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Abstract. 

We study the phase diagram of the extended Hubbard model in the atomic limit. 
At zero temperature, the phase diagram decomposes into six regions: three with 
homogeneous phases (characterized by particle densities p = 0, 1, and 2 and staggered 
charge density A = 0) and three with staggered phases (characterized by the densities 
p= ^,1, and ^ and staggered densities |A| = i, 1, and i). Here we use Pirogov-Sinai 
theory to analyze the details of the phase diagram of this model at low temperatures. 
In particular, we show that for any sufficiently low nonzero temperature the three 
staggered regions merge into one staggered region S, without any phase transitions 
(analytic free energy and staggered order parameter A) within S. 
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1. Introduction 

The theory of strongly correlated electron systems is nowadays a subject of 
vigorous research. The interest in these systems is stimulated, to a large extent, 
by attempts at explaining the mechanism of high-temperature superconductivity 
[MRR90, Dag94], the phenomenon of electron localization in narrow-band systems 
[IILM75] and properties of quasi one-dimensional conductors [Hub79], to name a 
few. Among the models that are most frequently studied is the Hubbard model 
augmented by a nearest neighbour interaction. This model, known as the extended 
Hubbard model, is defined by the following Hamiltonian 

Ht = -tY^ {cl^Cj,a + h.c.) + U ^ni,^ni,^ + W rurij - + zW +^)^ni. 

{i,j),a ieA ieA 

(1.1) 

In (1.1) we used the following notation: at each site i of a d-dimensional bipartite 
lattice A, with z nearest neighbours, there arc creation and annihilation operators 
c| ^ and Ci^cr of the electron with up and down spin, a I, that satisfy the 
canonical anticommutation relations, while rii^cr '■ = c| ^Ci^ a and ni : = rii^^ -\- rii^^. 
The first term of the Hamiltonian (1.1) stands for the isotropic nearest neighbour 
hopping of electrons, the second one is the familiar on-site Hubbard interaction, the 
third term represents the isotropic nearest neighbour interaction, and the last one 
the contribution of the particle reservoir characterized by the chemical potential 
We have introduced the shift zW -|- ^ in order to move the hole-particle symmetry 
point (the half- filled band) to the value /j, — 0. Originally, the second and the 
third terms were supposed to simulate the effect of the Coulomb repulsion between 
the electrons, hence only positive U and W were considered. Later on, in various 
applications of the model, the parameters t, U and W represented the effective 
interaction constants that take into account also other interactions (for instance 
with phonons). Therefore U and W could take negative values as well. In this 
paper U will be allowed to change its sign while W always stays positive. 

The so called narrow band case of the extended Hubbard model, i.e. \t\ ^ \U\, 
is of special interest in physical applications of the model. It has been studied by 
means of various approximate methods in many papers (see for instance [Lor82] 
and references quoted there). These studies revealed the existence of staggered 
charge-order at the hole-particle symmetry point. The staggered charge-order is 
characterized by a nonvanishing order parameter 

A= lim \A\-^yei{ni), (1.2) 

A/'oo ^ — ' 
ieA 

where £j assumes two values, 1 or —1, depending on to which sublattice of the 
bipartite lattice A the site i belongs, and (•) stands for the Gibbs state. 

Rigorously, the existence of staggered charge order has so far only been estab- 
lished in the so called atomic limit t — > [J§d94]. While the above mentioned 
approximate results suggest that the staggered charge order persists in the corre- 
sponding narrow band model, the methods used in [J§d94] unfortunately do not 
allow to establish this rigorously, because they rely on the reflection positivity of 
the atomic limit model which fails to be true for nonzero t. 

Here we propose to study the atomic limit of the model (1.1) using a different 
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also [Zah84,BI89]. On the one hand, these methods will allow us to study detailed 
properties of the low temperature phase diagram in the atomic limit, on the other 
hand they allow for an extension to nonzero t, treating the narrow band Hubbard 
model as a quantum perturbation of the t = model, and leading to the rigorous 
proof of staggered charge order in the narrow band Hubbard model [BK94]. 

In the atomic limit, it is convenient to rewrite the Hamiltonian (1.1) in a form 
that makes the hole-particle symmetry apparent. Namely, introducing Qi : = ni — l, 
we have 

{i,j) ieA ieA 

where we passed to dimensionless parameters U and fi, setting W = 1. Note that 
in the atomic limit all operators appearing in the Hamiltonian commute. There- 
fore, the model (1.3) can be viewed as a two-component classical lattice gas or, 
equivalently, as the classical gas with four possible states 0, t, i, 2 in each site, that 
correspond to an empty site, a singly occupied site with spin t or J,, and a doubly 
occupied site, respectively. 




FiG.l. Ground state phase diagram. There are three open regions 
Hq and with no staggered charge order (A = 0), and three 
open regions o}, (Si^ _}, and S*!- o}, with staggered charge order 
|A| = 1/2, |A| = i, and |A| = 1/2, respectively. 

In the sequel, we shall present and discuss the phase diagram of the model (1.3) 
on the lattice Z^. The ground state phase diagram is shown in Fig. 1. The (U^iJ.) 
plane decomposes into six open regions. In each of the regions H_, Hq, and H^, 
there is a unique homogeneous ground state, whose particle density p, given by 

^= A^J^r'E(^^+i)' (1-4) 
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equals 0, 1, and 2, respectively. The remaining three regions S'l.]. 'S'{_|_ o} and 
S'{_^o} are characterized by staggered charge-order. Namely, in the region S{a,b} 
there are two staggered ground states, denoted [a, b] and [6, o] , with Qi — a on one 
sublattice and Qi = b on the complementary sublattice. Note that the staggered 
order parameter A is nonvanishing in the whole staggered region S and jumps from 
|A| = 1 to |A| = I at the boundary between _} and 'S'j- o} or o}- 

Using reflection positivity it has been shown [J§d94] that the staggered long 
range order in the region S persists at small temperatures T = > 0, provided 
one stays sufficiently far away from the boundary between S and the homogeneous 
regions Ha, a — 0, ±. However, this does not exclude existence of a phase transition 
inside the staggered region S. Namely, in view of the discontinuity of A at the zero- 
temperature transition lines, one could expect that A reveals a phase transition at 
nonvanishing temperatures as well. Indeed, mean field arguments [MRC84] predict 
a first order transition surface emerging from the zero temperature transition line 
separating from •S'l+^o} and similarly for the line separating from 

Using a suitable notion of restricted ensembles we are able to analyze this ques- 
tion rigorously. Our main result here is to show the absence of any such phase 
transition, in contrast to the above mentioned mean field results. 

Theorem A. Consider the complement S of the union of closed homogeneous 
regions Ha, a = 0,±, 

S = R'^\{H-UHoUH+), 

and let 

= e ^ I dist(x, S^) > e}. 

For d > 2 and any e > there exists a constant (3o < oo (depending on e and d) 
such that, for all /Sq < /3 < oo and {U, n) e S^^\ there are exactly two phases^, 
(~)even ond (— )odd- Moreover, A > for the phase {—)even, A < for (— )oddj o.n'd 
both the free energy density, f{(3, U^y), and the staggered order parameters of the 
two phases, Aeyf;n{P,U, fx) and Ao^^{P,U, fx), are real analytic functions of U and 
H in S^"^^ . 

Rem,ark. As mentioned before, the zero temperature staggered order parameter A 
jumps from A = ±1 to A = ±^ at the boundary between S*!-!-^-} and 5'{_|_^o} oi' 
^i-^o}- It is interesting to relate this jump to the smooth behaviour for T = > 0. 
As we will see in Section 3, the crossover between these two behaviours is described 
by a hyperbolic tangent. Taking, e.g., the order parameter of the even phase in the 
vicinity of the boundary between, say, and 5'{+^o}) one obtains that 

3 1 , {2d- n-Ul2\ 

^ ~ T + T tanh — '— 

A A \ kT ) 

asT ^0. 

Turning to the homogeneous regions Ha, a = 0, ±, we use standard Pirogov-Sinai 
theory to discuss the low temperature behaviour inside the corresponding regions 
Ha^^ . This enables us to prove analyticity, unicity and translation invariance of the 
homogeneous phases. 



^As usual, a phase is defined as an extremal Gibbs state which is periodic in all d lattice 
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As for the boundaries between the staggered region S and the homogeneous 
regions Ha, we note that the zero temperature coexistence hne between Hq and 
gives rise to a coexistence hne surface of the two staggered phases with the 
homogeneous one. With decreasing (3, this surface bends towards negative U, i.e., 
above the ground state coexistence Hne between Hq and S^^^^y the corresponding 
homogeneous phase is stable [J§d94]. 

The remaining part of the zero temperature boundary between S and the homo- 
geneous regions is of similar type as the boundary between staggered and homoge- 
neous phases in the antiferromagnetic Ising model. For example, at the boundary 
between and H^, it is possible to join the staggered phase [+, — ], without 

any energy cost, to the second staggered phase [—,+], going through the homoge- 
neous phase stable in (see [J§d94] and also Section 2 below for a more detailed 
explanation). We therefore expect that this phase boundary turns, for nonzero 
temperatures, into a second order transition line. Actually, in the limit U ^ —oo 
the extended Hubbard model in the atomic limit becomes equivalent to the Ising 
antiferromagnet (with homogeneous external field fx) where this fact was rigorously 
proven [KY93] . In a similar way one expects that all other boundaries between ho- 
mogeneous and staggered phases, except for the boundary between Hq and 
already considered above, correspond to second order transition. 

We summarize our knowledge of the phase diagram of the extended Hubbard 
model in the atomic limit in the following theorem (see also Fig. 2). 

Theorem B. Let d > 2 and f3 be sufficiently large. Then there exist open regions 
Ha{(3), a = 0, ±, and S{(3), where Hq{P), and S{P) touch on a curve, 

L{l3)=dS{(3)ndHo{/3)j^(D, 

such that the following statements are true: 

i) For (t/, n) G S{/3) there exist exactly two phases, a phase (— )even with A > 
and a phase (— )odd with A < 0. These phases are periodic with period 2 
and both the free energy density, f{(3, U,ijl), and the staggered order param- 
eters of the two phases, Aeven(/?, t/, A*) cind Aodd{P,U, ji), are real analytic 
functions ofU and fx in S{l3). 

ii) For {U, n) G Ha{(3), a = 0, ±, there is exactly one phase (— )a- For this 
phase, A = 0, it is translation invariant, and the free energy density, 
f{(3, U, /J,), is a real analytic function of U and /x in in Ha{(3). 

iii) On the boundary L{f3) between S{j3) and Ho{P), the three phases {—)a, 
{~)even, o,nd (— )odd coexist. Furthermore, all periodic Gibbs states on this 
line are a convex combination of those three phases. 

iv) As ^^oo, HaiP) Ha, a = 0, ±, and S{P) S. 

Proof of Theorems A and B. Theorems A and B are immediate consequences of 
Propositions 1-4 that are stated and proved in Section 3. □ 

Before passing to the Propositions 1 - 4 we turn to a new representation of the 
model (1.3) in terms of spin 1 variables and then to a detailed examination of its 
ground state phase diagram. 
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/i/2c/ 




Fig. 2. Phase diagram at low temperatures. Thin lines denote 
the conjectured second order transitions, while the thick line, separating 
the two phase staggered region S from the homogeneous region Hq, is 
first order. Shaded are regions over which we have no rigorous control 
(they shrink to the zero temperature lines with /? ^ oo). 

2. Structure of ground states and restricted ensembles 

In order to rewrite the model (1.3) in terms of a classical spin system, we recall 
that all operators appearing in (1.3) commute. However, fixing all the eigenvalues 
Si of the operators Qi, Si e { — , 0, +}, does not completely specify the state of the 
system, because Si — corresponds to two possibilities n^^-f = 1 and rii^^ = or 
n^^i = and ^ — 1. In the partition function of the classical spin model, this 
leads to a factor of 2 for every singly occupied site, and therefore to an overall 
factor 2^*6A(i-<3i). in this way the quantum system (1.3) is mapped onto an 
antiferromagnetic spin 1 model, with 

H=J2SiS3 + ^T.^i-^'J2^i^ (2-1) 

where 

U = U-2p-^ln2. (2.2) 

It is useful to rewrite the Hamiltonian (2.1) as a sum over nearest neighbour terms 
h{Si,Sj), namely 

H=J2HSi,Sj\ (2.3) 
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where we introduced the energy per pair of nearest neighbour sites 

h{Si, Sj) = SiSj + ^{Sl + S]) -^{Si + Sj). (2.4) 

This form of the Hamiltonian makes the task of constructing the ground state 
phase diagram, i.e. determining the six regions Ha, a = 0, ±, S'{+^_}, S'{_|_^o}i and 
^'{-^o}) rnentioned in previous section, an easy exercise. The energies of the nearest 
neighbour spin configurations are 





2d d ' 


(2.5a) 


H+,0) = 


4- - ii 

^ ^ 4d 2d' 


(2.5b) 




-1 + 

^ 2d' 


(2.5c) 


MO,-) = 


r\ + K + JL 

" ^ 4d ^ 2d' 


(2.5d) 


h{0,0) = 


0. 


(2.5e) 


H-,-) = 


^ 2d ^ d ' 


(2.5f) 



Using (2.5), we find three homogeneous regions 

iy+ :={([/,//) I ^ > max{l, 1 + ^}}, (2.6a) 
with minimal energy pairs (+, +), 

H_ := {(t/,/x) I > max{l,l + ^}}, (2.6b) 
with minimal energy pairs (—,—), and 

Hq : = |(?7, ji) I max| 
with minimal energy pairs (0, 0); and three staggered regions 

with minimal energy pairs (+, — ) and (— , +), 

Sh,o,:={(C/..)I^<|^<1 + ^ a.d ^<A}, (2.7b) 
with minimal energy pairs (+, 0) and (0,+), and 

S(_,„,:={(C/,M)I^<-^<1 + |; and ^<-|j}, (2.7c) 

with minimal energy pairs (—,0) and (0,—). Thus, in each of the regions Ha, 



2d 



2}<4^l' 



(2.6c) 



2d 
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other hand, in each of the regions S^a,b} there are exactly two ground configurations, 
such that, when on one sublattice Si = a, on the complementary one Si = b, and 
vice versa. In order to relate the formulae (2.6) and (2.7) to Fig. 1, we notice that 

f/ = C/ for /3 = oo. 

At this moment, let us remark that the above analysis of ground configurations 
shows that the components h of H (cf. (2.4)) constitute a so called m-potential 
[Sla87]. Moreover, since in each of the regions Hq, H_, /S'{+^_}, -S'{+^o}) ^^nd 
(Sj- 0} there are only finitely many ground configurations, one can readily apply 
standard Pirogov-Sinai theory to study the low temperature properties of the cor- 
responding phases, away from the boundaries of these regions. 

On the boundaries of the regions Hq, S^^^^y, o}, and o}, the 
structure of the ground states is more complicated. Combining the minimal en- 
ergy pair configurations of the corresponding adjacent regions, one can construct 
infinitely many ground configurations everywhere, except for the boundary between 
S'^-i-^-}. and Hq. On the latter boundary there are exactly three ground configura- 
tions, namely, those that correspond to the adjacent regions. This, of course, places 
also this case into the realm of standard Pirogov-Sinai theory. 

There is an important difference between the infinitely degenerated boundaries 
shared by staggered and homogeneous regions and those shared by two staggered 
regions. In the first case, i.e. on the boundary between, say, and S*!-!- the 
minimal energy pairs of both regions, namely the pairs (-I-, -|-), (+, — ) and (—,+), 
can be combined into an arbitrary configuration made out of "-h" and "— ", as 
long as no nearest neighbour pair of minuses is present. Mutatis mutandis, the 
same is true for the other infinitely degenerate boundaries between staggered and 
homogeneous regions. One therefore obtains the same structure of ground states as 
in the Ising antiferromagnet at the critical field, presumably giving rise to surfaces 
of second order transitions at nonzero temperature. 

In the second case, i.e. on the boundaries shared by two staggered regions, 
the situation is different. Considering for instance the region S+ that consists 
of S"!-!-^-}, S'{_|_^o}: and the boundary shared by these regions, we introduce two 
disjoint classes of configurations, ^J^gn ^^^^j. Namely, we define Q^y^n ^ ^^e 
set of all configurations, for which 5'i = -|- on the even sublattice, while Si = 
or — on the odd sublattice, and the set G^^^ by interchanging the role of the two 
sublattices. Then, for all points in S+, every groimd configuration falls into one of 
these two disjoint classes. This already suggests to use, in 5-1-, a version of Pirogov- 
Sinai theory with the sets ^J^gn ^^^^ ^odd P^^yi^S the role of restricted ensembles 
[BKL85] that replace the ground states^. The same remark applies to the region 
5'- made of S'l^. 'S'{_ o} and their shared boundary, and the corresponding sets 

Geven and G'^^. 



^Notice that different configurations from ^^en (and similarly G^^^) attain, in general, different 
energies. In particular, only on the boundary between S'{_|_ _j and q}, all configurations from 
Seven are ground configurations. As a consequence, the version of the Pirogov-Sinai theory for 
lattice systems with residual entropy [GS88, SGL89] does not apply here. Namely, it would need 
in our case an assumption that everywhere in the coexistence region S+ all configurations from 
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3. Proof of Theorems A and B 

As mentioned above, Theorems A and B are an immediate consequence of the 
following four Propositions. 

Proposition 1. Consider the regions 

= {xeHa\ dist(a;, Hi) > e}, a = 0, ±. 

For d > 2 there exists a constant b = b{d) > 0, such that for all e < oo, j3 > h/e, 
and {U, e Ha \ there is exactly one phase {—)a- For this phase A = 0, it is 
translation invariant, and the free energy density /(/?, U, fj,) is analytic in 

«<.) = {(p,U„, e C3 1 R« , > V. (^If . ^) e ^<'.}. (3.1) 

Proof. Except for the last statement, the proposition follows immediately from 
standard Pirogov-Sinai theory [PS75, Mar 75, Zah84]: Given the relations (2.5) 
and (2.6), one gets, for a suitable constant a = a{d) > and {U, fi) e Hi^\ the 
inequality 

Ph{b,c) > Ph{a,a) + f3ae (3.2) 

provided (&, c) ^ (a, a). As a consequence, all excitations of the ground state (a, a) 
are exponentially suppressed, with a "Peierls constant" r = d^e where o: > 
depends only on the dimension d. 

In order to prove the last statement of the proposition, one needs a representation 
in terms of a convergent contour expansion inside TC'^\ where the Hamiltonian H 
is complex. This situation has been dealt with in [BI89]. In order to prove the 
corresponding Peierls condition, one needs a relation of the form (3.2) for the real 
part of (3H, namely 

Re ph{b, c) > Re ph{a, a) + (Re p)ae . (3.3) 
This leads to the regions H^'a^ . □ 

Proposition 2. Consider the regions S!^a,b} o,i^d the corresponding sets 
{a, 6} = {+,—}, {+,0}, { — ,0}. For d> 2 there exists a constant b = b{d) > 0, 
such that for all e < oo, P > b/e, and {U, fi) G b}' ^^^^^ exist exactly two phases, 
the phase (— )even with A > and the phase (— )odd with A < 0. These phases are 
periodic with period 2 and the free energy density, f{(3, U, (i), as well as the corre- 
sponding staggered order parameters, Aeven(/5, U, //) and Aodd(/5) U, jj), are analytic 
in 

= {(A e I , , V. (^, If) e .it,}. (3.4) 

Proof. Again, the proof is standard, except for the last statement. Actually, in view 
of the essential singularities associated with first-order phase transitions [Isa84] , the 
analyticity proof in the coexistence regions S^^ab} ^oie subtle than that in the 

„; ,!„ „1 
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We start with the observation that a Peierls condition of the form (3.3), namely 
Re ph{c, d) > Re ph{a, b) + (Re P)ae for all {c, d} ^ {a, b} , (3.5) 

is valid in all of S^^^^y. Introducing truncated contour models as in [Zah84, BI89] to 
expand about the two staggered ground states [a, b] and [6, a] , one therefore obtains 
a convergent cluster expansion for the corresponding "truncated free energies" /even 
and /odd- Next we note that 

/even(/?, U, fi) = fodd{P, U, fj.) for aU (/?, U, fi) e S['l,y 

due to the translation symmetry of the model. As a consequence, 

Re (/3/even(/3, U, fj)) = Re (/3/odd(/3, U, fj)) for aU (/3, [/, //) e ^{^j,^ . 

The results of [BI89] then imply that both the even and the odd phase are stable 
in all of S^^^ , implying in particular that the truncated free energies are equal to 
the corresponding "physical free energy" /(/?, [/, //) obtained as the limit of (log- 
arithms of) finite volume partition functions. As a consequence, /(/3, U, //) can 
be expressed as an absolutely convergent sum of analytic terms, implying that 
/(/3, U, n) is analytic itself. The stability of both the even and the odd phase also 
implies the convergence of the contour expansion for the staggered order parameters 
Aeven(/5, U, jj) and Aodd(/5, U, ji), yielding their analyticity in <S|^^ j,}. □ 

Remark. Let us note the differences between the situation of Proposition 2 and the 
phase coexistence of, say, an Ising ferromagnet at zero field h. In the situtation 
of Proposition 2, where the two phases (— )even and (— )odd can be obtained from 
each other by a translation, /even = /odd troughout the complex region S^^b}^ ^ 
statement which is stronger than the stability condition Re {(3 f even) = Re (/3 /odd)- 
For the Ising model, on the other hand, the symmetry relating the two phases 
(— )+ and (— )- requires a change of the sigh of h. As a consequence, no open 
neighbourhood U C C oi h = can be found such that both the plus and the minus 
phase are stable in U. Furthermore, on the coexistence line Re h = where both 
phases are stable, f+ ^ f^, even though the weaker condition Re (/?/+) = Re (/?/-) 
is true for Re h = 0. 

Proposition 3. Let R be the union of Hq, _} and their common boundary, and 
let R^^^ = {x E R \ dist(a;, R'^) > e}. For d > 2 there exists a constant b = b{d) > 0, 
such that for all e < oo and (3 > b/e there exist a curve L{(3), separating R^^^ into 

two open regions R^^^ and _y such that the following statements are true. 

i) In R^q \ there exits exactly one phase (— )o- This phase is translation invariant 
with A = 0. 

ii) In _| there are exactly two phases (— )even and (— )oddj characterized by 
A > and A < 0. Both phases are periodic with period 2. 

Hi) On the boundary L{f3) between R^^ and i?^^ all three phases (— )even 

and (— )odd coexist. Furthermore, all periodic Gibbs states on this curve are a 
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iv) As a function ofU and fjb, the free energy f is real analytic in R^'^^\L{f3), and the 
staggered order parameters of the two phases (— )even and (— )oddj ^eveniP,U, /i) 
and Aodd(/5) U,ii), are real analytic in i?^^ _y. 

Proof. Again the proof is standard. One now introduces three different truncated 
contour models: one for the excitations about the homogeneous configuration (0, 0) 
and two for the excitations about the two staggered configurations (+, — ) and 
(— , +). In the region R^'^\ and more generaUy in the complex region 

7^(^) = {{U,ii) I (Ret/, Re//) G R^'^ , (3.6) 

the model again satisfies a suitable Peierls condition provided /?e is big enough. This 
leads to the convergence of the cluster expansion for the corresponding truncated 
free energies /o, /even and /odd in 7^(^) D R'^^\ 
Given the "degeneracy removing condition" 

A(,(+,-)-MO.O))=i>0 (3.7) 



and the symmetry relation 



/even([/, fl) = /odd(t^, 1^) , (3.8) 



the proof of statement i) - iii) is now an easy application of the methods of [Zah84] . 
Actually, the complex analogue of (3.7), namely the degeneracy removing condition 



d 



dReU 



(^Re/i(+,-) -Re/i(0,0)^ = ^ > 0, (3.9) 



together with the validity of (3.8) in the complex region TZ^^'' leads to the existence 
of a phase transition surface S{P) that separates TZ^^^ into two open regions: a 
region TZ'^q^ in which Re{Pfo{U,n)) < Re (/3/even(t^, /^)) and f{U,ii) = fo{U,n), 
and a region TZ^^^ _y in which Re (/3/o(C/, //)) > Re (/3/even(t^, A*)) ^^nd f{U,ii) = 
feven{U, fi), See [BI89]"^. As a consequence, the free energy / of the model can be 
rewritten as a convergent sum of analytic terms in both T^.q^'* and TZ^^^ _t., leading to 

the analyticity of / in TZ^^^ \ jC{/3) and hence the real analyticity of / in R^^^ \L{/3). 
In a similar way, one obtains the real analyticity of the charged order parameters 

^even 

Proposition 4. Consider the regions S± introduced in the last section, and the 
corresponding sets S^^ . For d > 2 and m = ± there exists a constant b = b{d) > 
0, such that for all e < oo, f3 > b/e, and {U, n) e 5'^'*, there exist exactly two 
phases, a phase (— )even with A > and a phase (— )odd with A < 0. These 
phases are periodic with period 2 and both the free energy density, /(/?, U,fi), and 



^In the lanffuaffe of [Zah84, BI89], Tz''^'^ is the region where the homogeneous phase is stable, 
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the corresponding staggered order parameters, Aeven(/3, U, /i) and Ao^^{l3, U, /i), are 
real analytic functions of U and jJL in Srn . 

Proof of Proposition 4- Without loss of generality, we may assume that {U, /j,) e 
S+. In order to prove the proposition, we introduce an auxiliary Ising variable 
cTj = (7{Si) as 

and rewrite the model (2.1) in terms of an effective Hamiltonian H^^{a). We then 
prove that the Hamiltonian H^^ has two ground states g^"^^"^ and (7°'^'^, correspond- 
ing to the restricted ensembles ^even ^odd introduced at the end of the last 
section, 

^r" = ei and = -e^ , (3.11) 

and that the excitation above these ground states obey a suitable Peierls condition. 
Here, as in Section 1, = +1 on the even and ti = —1 on the odd sublattice. 

We start with some notation. We consider a box A = [1, L]^ n Z^, its boundary 
dh. = {z e A*^ I dist(z. A) = 1}, the set 5(A) of nearest neighbour bonds (i, j) with 
at least one endpoint in A, and the union of A and its boundary, A = A U dK. Here, 
as in the sequel, dist(-, •) denotes the £1 distance in Z'^. As usual, we call two sets 
F, V' e Z*^ adjacent or touching if dist(V, V') = 1, and a set F C connected if 
for any two points i,jEV there is a sequence of adjacent points in V that joins i 
to j. 

Keeping in mind that we want to construct finite temperature states (— )m which 
are small perturbations of the restricted ensembles m — "even" or "odd", we 
introduce an effective Hamiltonian H^{aA | m) in A with the boundary conditions 
m = even, odd, by 



cT(Si)=ai,ieA <i,j)eB(A) 



Sa- 

The corresponding finite volume Gibbs states are 

(■>™.A=Z^E--^''^"<-'"" (3-13) 

^ ^ CTA 

with 

CTA 

Extending the configuration ua to A by setting ai = g^^ for i e dA, we define a 
nearest neighbour pair G B{A) as excited in the configuration a a if ai = aj 

and a point i G A as excited if it is contained in an excited bond. Note that the 
notion of whether a bond that joins the volume A to its boundary dA is 

excited or not depends on the boundary condition. 

At this point, contours and ground state regions are defined in the standard way: 
Given a configurations a\ (and one of the two boundary conditions introduced 
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as pairs of the form Y = (suppF, ay), where suppF is a connected component of 
the set of excited points and cry is the restriction of cr^ to suppY. The ground 
state regions are defined as the connected components of the set of points which are 
not exited. Note that the restriction of to a ground state region C is staggered, 
and hence equal to the restriction of one of the two ground states to C, ac — gci 
where m = m{C) may vary from component to component. 

An important property of a set of contours li, • • • corresponding to a con- 
figuration (Ta is that they "match". In order to define this notion, we note that 
each contour Y determines the value of cr^ on all points i G A which touch its 
support because all bonds joining the support of Y to such a point are not exited. 
The contour Y therefore determines the value of m(C) for all ground state regions 
C touching its support. We say that Y attaches a label m{C) = my(C) to these 
ground state regions. Matching of the contours li, ■ ■ ■ jY^ is the statement that the 
labels attached to a given ground state region C by different contours are identical 
and compatible with the boundary condition. A minute of reflection now shows 
that to each set {Fi, • • • , F„} of matching contours with dist(supp Y^, supp Y/) > 1, 
k ^ I, there corresponds exactly one configuration a\. The partition function 
Z^(A) can therefore be expressed as a sum over sets of matching contours, once 
the Hamiltonian H'^{a\ \ m) has been expressed in terms of Yi, • • • , Y^. 

We will now show that this can be done in the form 



n z{Yu) (3.15) 



g-^ifr(aAlm) ^ g-^Hr(3r|m) 

where ziY]^) are contour weights obeying a Peierls condition 

< e-"l'"PP^I (3.16) 

with sufficiently large Peierls constant r. 

We start with an explicit calculation of the Hamiltonian HJ^{a\ \ m) for the 
configuration cta = g]^ with no contour. In this configuration, each point i e A 
with cTj = — has 2d nearest neighbours j E A with ctj = + if « e A, and 1 nearest 
neighbour j & A with aj = + if i E dA. Since aj = (T{Sj) = + implies Sj = +, the 
summation over the spin variable Si in (3.12) therefore leads to a factor 

X= ^ g-2d/3h(Si,+) ^ g-2d/3/i(0,+) _^ ^-2df3h{-,+) ^3 -^7^ 

Si:a{Si)=- 

if i e A, and to a factor 

Si:a{Si) = - 

if i e 9A. Per bond, this yields the energy 

r h'o = -4 log A' ifiedA and ai = g^" = - 
/i-((i,j)) = J . (3.19) 

[ ho = -^\ogX if i e A and ai = gV^ = - 

For the energy of the configuration g^, this gives 

Hf{gT\m)= h^i^^^j)). (3.20) 
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Remark. Obviously, the boundary correction (3.19) does not affect the specific 
ground state energy 

em = lim T^MfigT I m) (3.21) 

A^Z<i |A| 

so that Ceven = Codd- It does affect, however, the finite volume ground state energies 
HfigT I m). 

In order to calculate the weight e~^^^ (o-a|"i) ^^j, configuration a a correspond- 
ing to a nonempty set of contours {Yi • • • , F^}, we extract from (3.12), for each 
contour Y e {Yi ■ • ■ ,Yn}, the factor 

z{Y)= n e-f'^^^^'^^K (3.22) 

SY:aiSi)=ai {i,j)eB{Y) 

Here B{Y) is the set of all bonds G -B(A) such that either 

i) both endpoints of are in the support of Y, 

or 

ii) only one endpoint of {i, j) lies in the support of F, and this endpoint corresponds 
to a value cr^ = —1. 

Note that the second class of bonds are those bonds which couple the spin variables 
in the support of Y to the spin variables in A \ supp Y . The remaining sum in 
(3.12) can be easily calculated because all points i G A \ (supp Yi U ■ ■ ■ U supp y^) 
with (Ji — — are not excited. The summation over the corresponding spin variable 
Si therefore again leads to factors A and A', giving a factor e-^'»"«*'J)) for all the 
bonds in B{K) \ {B{Yx) U • • • U B{Yn)). Extracting the factor 

-Q g-/3ft^(({ij))) 
{i,3)eB{Y) 

from the activities (3.22), we therefore obtain a representation of the form (3.15), 
with 

z{Y)= n e-^(^(^-^^)-^o-((iJ))) . (3.23) 

SY:^{S^)=ai {i,j)eB{Y) 



We are left with the proof of the Peierls condition (3.16). We start with the 

+ 



observation that for (U, jj) e C 5"+ and arbitrary values for the spins Si and 



h{Si, Sj) > min{/i(0, +), h{-, +)} , (3.24) 

while 

h{Si, Sj) > min{/i(0, +), h{-, +)} + ae , (3.25) 

for some dimension dependent constant o: > whenever the bond is excited. 

Combining (3.24) and (3.25) with the fact that 

/i^((z, j)) < min{/i(0, +), h{- +)} , (3.26) 

we obtain the bound 

\<Y)\< Y n (3-27) 
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where B*{Y) is the set of excited bonds in B{Y). Bounding from below by 

^1 suppy|, and the number of terms in the sum over Sy from above by 2I^"pp^I, 
we obtain the bound (3.16) with 

e-^ = 2e-5/3'*^ (3.28) 

Given (3.14), (3.15) and (3.16), the partition function (3.14) can be expressed 
as the partition function of a contour system with exponentially decaying weights, 

n 

Z^(A) = e-/3^A"(^A l™) J2 U'(^^) (3-29) 

{Yi,---,y„} k=i 

where the sum runs over sets {Yi, ■ ■ ■ ,Yn} of matching contours obeying the com- 
patibility condition that dist(supp Yfe, supp Y/) > 1 for /c 7^ L As a consequence, the 
model can be again analyzed by standard methods, see e.g. [Zah84]. One obtains 
the existence of the limits 

{-)m = lim (-)^,A (3.30) 

L— »oo 

as periodic Gibbs states with A > and A < 0, respectively, the fact that these 
states are extremal, and the fact that all other periodic Gibbs states are convex 
combinations of (— )even and (— )odd- 

Considering finally a suitable complex neighbourhood of the region S^^ , e.g. 

S^^'^^ = {(U, 11) I (Re U, Re /i) e 5^'^ and |Im t/| < S, |Im ii\ < 5} (3.31) 

with 5 sufficiently small, one easily establishes a Peierls condition with a slightly 
smaller Peierls constant r = r(/3, e, 5, rf). The methods of [BI89] then give'' the 
free energy density / and the staggered charge-order parameters Aeven(/5, t^, At) 
and Aodd(/3, t^, /u) as convergent sums of analytic terms in S^'^\ implying their 

analyticity in S^'^^ and hence their real analyticity in S^K This completes the 
proof of Proposition 4. □ 

Remark. It is intriguing to relate the first order jump of the staggered order pa- 
rameter at T = to the analytic behaviour at positive temperatures. To this end, 
we note that the distribution of the spin variable Si in the restricted ensembles 
is given by 

if 9? = + 

[- A {S{Si,0) + 6{Si,-)) if gT = - 

The finite temperature excitation above the corresponding ground states a = g'^ 
slightly modify these distributions, leading to corrections of the order 0{e~^'^^). 
As a consequence, the staggered order parameter A in the region S^^ is given by 

A = ±(^^ + ^ tanh(^2d/3(/i(0, +) - h{-, +)))^ + ©(e"^"^) 

= ±{j + l tanh(/3(2d - - C//2))) + 0(e-'3"^) , (3.32) 

where the plus sign corresponds to the even phase and the minus sign corresponds 
to the odd phase. 



■^As in the proof of Proposition 2, we use the fact that, by translation invariance, the corre- 

• J. _1 J" 1 • (->(C,<5) 
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